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ABSTRACT 



The vestiges of planet formation have been observed in debris disks harboring 
young and massive gaseous giants. The process of giant planet formation is 
terminated by the dissipation of gas in the protoplanetary disk. The gas-rich disk 
around HD 142527 features a small inner disk, a large gap from ~10 to ~140 AU, 
and a massive outer disk extending out to ~300 AU. The gap could have been 
carved-out by a giant planet. We have imaged the outer regions of this gap 
using the adaptive-optics camera NICI on Gemini South. Our images reveal that 
the disk is dynamically perturbed. The outer boundary of the roughly elliptical 
gap appears to be composed of several segments of spiral arms. The stellar 
position is offset by 0.17 ± 0.02" from the centroid of the cavity, consistent with 
earlier imaging at coarser resolutions. These transient morphological features 
are expected in the context of disk evolution in the presence of a perturbing 
body located inside the cavity. We perform hydro-dynamical simulations of the 
dynamical clearing of a gap in a disk. A 10 Mj up body in a circular orbit at 
r = 90 AU, perturbs the whole disks, even after thousands of orbits. By then 
the model disk has an eccentric and irregular cavity, flanked by tightly wound 
spiral arms, but it is still evolving far from steady state. A particular transient 
configuration that is a qualitative match to HD 142527 is seen at 1.7 Myr. 

Subject headings: Protoplanetary disks — Planet-disk interactions — Stars: individual 
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Introduction 



The lifetime of disks sets the time available for the planet formation process. In the 
Solar System, chronology of the oldest solid components based on radionuclides (e.g., Al 26 ) 
i ndicates that the dura tion timescale of the Solar Nebula was short, of order 2-3 M yr 



(IMontmerle et al. 



20061 ) . Direct imaging of a 9 ± 3 Mj up planet (ILagrange et al. 



20_nj) inside 



the debris disk of j3 Pic proves tha t large gas giants can form by I2I4 Myr. In protostellar 



systems, by 5-10 Myr no gas is left (IPascucci et al. 



20061 ) . The evidence points at early gia nt 



plane t formation, as in the candidate ma ssive protoplanets LkCa 15 b (IKraus &: Ireland 



2011) (although both remain to be confirmed, e.g. by 



20121 ) and T Cha b flHuelamo et al. 
direct imaging). 

Radial gaps in gas-rich protoplanetary disks are currently thought of as possible 



signposts of planet formation. Previous H and K 



)and coronograp 



reveal a hole in the disk, about 100 AU in radius (IFukagawa et al. 



lie im ages of HD 142527 



is in fact a gap with an inner boundary, abutting on an inner disk ( 



2006 1). This inner cavity 



van Boekel et al 



rough ly 10 AU in radius (Anthonioz et al., in preparation), and u p to 30 AU (|Verhoeff et al. 



2004), 



20111 ). that accounts for the large near-IR excess of HD 142527 (IFukagawa et al 



2010). 



Radiative transfer modelling of the spectral-energy distribution ( SEP) and NIR imag es are 



consistent with an inclination angle of ~20 deg (close to face-on, 



Verhoeff et al. 



201 lh . 



The youth of HP 142527A, at ~ 140 pq£|, is evid ent from the copiou s amounts of gas 



Oberg et al 



2011 



) in its ~0.1 M circumstellar disk (jVerhoeff et al. 
(120061 ) estimate a stellar age of 2^ 2 Myr and a mass of 1.9±0.3 M G 



2011) 



Fukagawa et al. 



©■ 



Here we report on deep adaptive optics near-IR imaging of the disk and gap of 



1 Assuming it is associated to Sco OB2; the Hipparcos parallax of 4.29±0.98 marcsec 
lower-limits its distance to 138 pc, at 3 o. 
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HD 142527. § |2] describes our observational setups, § [3] emphasizes our main observational 
results, § S] interprets, and § [5] concludes. 



Observations 



NICI adaptive optics imaging. We conducted 



HD 142527 with the NICI coronograph (IBiller et al. 



R ad aptive optics observations of 



2010) on the Gemini telescope. NICI is 



a dual-channel imager, the images were acquired simultaneously, by pairs, a different filter 
in each channel. Here we report on two pairs of filters: 1- broadband Ks (2.15 /an) with 
H 2 (1-0)S(1) (2.12 //m), both with a 0.22" mask occulting 95% of the stellar light, and both 
exposed to 2 s and 30 s, for a total of 1 h; 2- broadband L' (3.8 /im) with Br7 (2.17 /im), 
both maskless and exposed to 0.38s, for a total of 1 h. All setups were acquired in stare 
mode, compensating for parallactic angle rotation. The result of our NICI imaging is shown 
in Fig. [U after PSF subtraction^. 

SINFONI integral field spectra. We also conducted complementary integral-field-unit 
spectroscopy with the SINFONI instrument on the VLT in order to search for possible 
planetary objects in the gap. In our setup we chose a very narrow jittering throw that 
samples the inner 0.45 x 0.45 arcse c" 2 only. We proces sed the ESO pipeline datacubes using 



spectral deconvolution techniques (IThatte et al 



20071 ). and also through conventional PSF 



subtraction (using a P SF standard s t ar). The detection limits we achieve are consistent 



with those achieved by 



Thatte et al. 



( 120071 ) in AB Dor with the same instrument (~9 mag 



contrast at 0.2 " separation). 



2 we used a Moffat profile fit to the wings of the PSF 
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Observational results 



3.1. Size and morphology of the cavity and outer disk inner rim 



The NICI images clearly reveal the non axi-symmetric shape of the cavity. This is 
the highlight of our observations. The shape is roughly elliptical, albeit with significant 



departures, with a representative eccentricity of e c = yl — (b c /a c ) 2 ~ 0.5 (where a c and 
b c are the ellipse axes). But, as illustrated in Fig. [2j the exterior ring can be described 
as being composed of four arms overlapping in their extremities, except for two small, 
~0.2 arcsec gaps or intensity nulls along the ring, due north and south-south-east from the 
star. The observation of these morphological feature is based on the shape of the inner rim 
of the ring seen in direct images (even without subtraction of the PSF glare). Additionally, 
an unresolved brightness increment, or knot, can be seen along the north-western segment 
(which is also the brightest segment) - the alignment of this unresolved feature in all filters 
is illustrated in the RGB image, although we cannot determine from these data alone if this 
knot is real or perhaps only a ghost or a chance alignment of speckles. 

The south-western spiral arm segment (labelled arm 2 on Fig. [2]) sprouts away from the 
outer ring at a PA of - 45 deg (East of North), where it extends into the outer arm seen by 



Fukagawa et al. 



( 120061 ) . This outer arm is at much fainter levels than the outer ring - our 



dataset has finer resolution but is s hallower, so that only its root is seen in our figures. As 



discussed in 



Fukagawa et al. 



(120061 ) . this outer arm could have been triggered by a stellar 



encounter, but the putative p artner remains to be i dentified. A nearby point source turned 



out to be a background star (jFukagawa et al 



20101 ) - in our NICI images this source is 



located at a separation of 5. "44 from HD 142527, and a PA of 220 deg (East of North). 

The inner radius of the ring varies from 0.7" to the north, to about 1.0" to the south, 
and a projection effect is discarded since the central star is clearly offset by 0.17 ± 0.02" 
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from 



;he centre of this approximate ellipse. This offset was first noted by iFukagawa et al 
( 120061 ) . albeit at lower angular resolution. 



In if-band the w e st side of the rim we see is brighter, in agree ment with 



Fukagawa et al. 



(120061 ) . 



Fujiwara et al. 



(I2006h and 



Verhoeff et al. 



(120111 ) already noted that 



in the thermal IR, TV an d Q bands, it is the e astern side that is brightest, which is consistent 



with a projection effect (IFujiwara et al. 



20061 ). At the shorter wavelength i^-band, the front 



(i.e. nearest) side of the disk appears brighter because of forward scattering at the surface 
of the disk. On the opposite, the thermal emission from the larger projected surface of the 
rim of the back side of the disk appears brighter than the slimmer (in projection) rim of the 
near side. This implies a major axis of the disk in the North-South direction and a minor 
axis in the E-W direction (as observed here). The emission null we detect due North is seen 
on both sets of Subaru images, at f^-band and in the mid-IR as well. 



3.2. Detection limits in the range 14-35 AU 



No companion is clearly detected in our NICI images. An upper limit on a binary 
companion is difficult to obtain because of the speckle noise, and because of the lack of 
accurate NICI zero points. However, the SINFONI observations allow us to place a lower 
limit flux ratio at 0.2-0.3" separatio n of 1905 at 3 a in K; Fig. [3] illustrates this upper limit. 
Modeling the SED (with MCFOST, 



Pinte et al. 



20061 ) we estimate that half of the total 



flux density of HD 142527 in K can be assigned to the disk. In that case, the SINFONI flux 
ratio limit of 1905 corresponds to an absolute magnitude Mk > 8.2 at 3 a. As a comparison 



point we refer to the C0ND03 tracks ( iBaraffe et al 



20031 ) at 1 Myr of age, which would 



imply a mass of less than 12 Mj up given this limit. However the age of HD 142527 and 
the models are not precise enough to ascertain this upper limit mass value. We stress again 
that this limit constrains the presence of massive gaseous giants at stellocentric radii of 14 
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to 35 AU. 



4. Analysis 



4.1. Arguments in favor of a planet-created gap for HD 142527 



The contribution from photoevaporation to cavity clearing can be neglected in this 



case. The inner disk, the large gap and the vigorous accretion rate (IGarcia Lopez et al. 



20061 ) of 7 x 10 -8 M yr _1 argue against photoevaporation. An X-ray driven wind can 
indeed explain a subset of the observed accreting disks with inner holes , but only those w ith 
radii smaller than 20 AU and accretion rates less than 10 -8 M yr _1 ( jOwen et al.ll201l[ ). 



Gap clearing by single or multipl e planets could potentially explain wide gaps observed 



in transitional disks ( Mice et al 



2003|). We thus focus on the simplest possible scenario to 



interpret the morphology of the cavity in HD 142527's disk, i.e. dynamical clearing by one 
single embedded and massive- Jupiter-size planet opening a gap in its gaseous disk. The 
maxim um gap opening width will depen d on the mass and the orbital parameters of this 



planet (IDodson-Robinson fc Salvk 



the star-planet binary (lArtymowicz fc Lubow 



2011 



). The 



cavit y forms due to the torque exerted by 



19941 ) . Elliptical cavities have been predicted 
with e c ~ 0.25 in the case of circular binaries with mass ratios > 3 x 10 -3 , in which case 



the gap becomes large enough to deplete the region in which t 



outer 1:3 Lindblad reso nance is produced ( IKley fc Dirksen 



( iHosseinbor et al. 



le eccentricity-damping 



20061 ). Further numerical work 



20071 ) shows that, for eccentric binaries, the outer edge of the cavity 



becomes even more eccentric. The observed eccentric cavity is then a recurrent feature of 
gap- clearing models. 
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4.2. Fargo hydrodynamic simulations of gap clearing by a planet 



In order to test further the single-planet origin for the wide gap in the gaseous disk of 
HD 142527, we conducted hydrodynamic simulations with FARGO. Our goal here is not 
to fit the data but to verify qualitatively the tenability of this scenario fo r HD 142527 . 



20001) 



FARGO is a dedicated grid-based 2D code publicly available on the web (IMasset 
and specifically designed for planet-disk interactions. The run presented in this paper 
corresponds to a system evolution equivalent to 1.7 Myr. 

We picked common fiducial parameters for the model, as detailed below. The 
model disk is initially axisymmetric, its surface density is initially a power law of radius, 
£ = £ (r/ro)~\ where E = 566g cm" 2 (6.3 x 10~ 5 AIT 2 ) at r n = 1 AU, so slightly 



smaller than the value assumed by 



Dodson-Robinson fc Salykl (120111 ). The model disk has 



a constant pressure scale height H over radius r with value H/r = 0.05. For the disk 
calculations the boundaries are 4 and 300 AU for the inner and outer radii, respectively, 
giving a total mass of 0.1 M Q . The inner boundary is open, allowing mass to fall onto the 
star, while the outer boundary is such that we allow for a steady state mass transfer into the 



disk. Accretion across the disk is modele d by using an a prescription (IShakura fc Sunyaev 



19731). We adopted a = 0.002 following 



Dodson-Robinson fc Salykl (120111 ). From the 



SINFONI limits described above we know that any planet present in the inner 14 to 35 AU 
cannot have a mass exceeding ~12 Jupiter masses. We assume no such planet is present 
further out in the disk either, and chose a planet mass of 10 Mj up . 

We placed a 10 Mj up planet in a fixed circular orbit at 90 AU in the model disk. This 



radi us was chosen so that the p 



le.g 



Dodson-Robinson fc Salyk 



anet's sphere of influence, given by a few times its Hill radius 



20111 ). reaches the observed gap size of ~130 AU. Migration 



was inhibited until a sizeable gap was carved by the plan et. This artificial constraint stem s 



from the lack of heat or momentum diffusion in FARGO (jDodson-Robinson fc Salvia 



201lh . 
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The resulting gas surface density field is shown in Fig. |U This specific dynamically 
perturbed morphology will persist over a timescale of a few tens of thousand years (i.e, 
around 100 orbits) but is not in steady state. The rim will remain asymmetric over 
time, and is composed of tightly wound spiral arms, whose superposition forms a roughly 
elliptical cavity. The cavity centroid is persistently offset from the star, but any particular 
morphology of the gap is transient and lasts about 100 orbits. In this snapshot view (fig. HJ), 
the outer edges of the cavity are most reminiscent of those seen in the NICI images. The 
cavity eccentricity is ~0.6, and its centroid is offset from the star by ~16 AU. 

As seen in Fig. m this model produces a massive inner disc, extending out to 50-60 AU 
in radius, which is not consistent with the observations. Additional ingredients are required 



to deplete that region in the fr amework of a dynamical c 



planets closer to the star (as in 



Dodson-Robinson fc Salyk 



earing model, e.g., additional 



20111 ). Fitting the inner disc is 



left for future work. We stress that the purpose of these simulations is to provide a plausible 
explanation of the NICI observations of the inner edge of the outer disk only. 



5. Summary and discussion 

Our observations of HD 142527 indicate that the inner rim of its outer disk is roughly 
elliptical, and composed of an overlapping set of spiral arm segments. We confirm that the 
star is offset by 0.17 ± 0.02" from the centroid of the cavity. No companion is detected in 
the range of 14 to 35 AU, with a lower limit Mx > 8.2 at 3 a (for a distance of 140 pc). 

The observed morphology of the disk in HD 142527 suggests a dynamically perturbed 
state. The segmented spirals, and the offset between the cavity centroid and the star, 
can be seen among the varied morphologies predicted by hydrodynamic simulations of 
HD 142527. To explore, at least qualitatively, the origin of the cavity and its peculiar shape, 
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we reported on a simple case. A massive protoplanet on a circular orbit whose radius is 
comparable to the disk's outer edge perturbs the entire disk, which does not reach steady 
state even after thousands of orbits. The global shape of the cavity is a long-standing 
feature; the general asymmetries (off-center star, eccentricity of the cavity) are imprinted 
early and remain for the duration of the calculations. However, the specific configuration of 
the model run, position of the streamers, and position and number of spiral arms, form a 
specific morphology that is only a transient phase. A particular morphology may recur, but 
lasts for a limited number of orbits. That phase which qualitatively matches HD 142527 
lasts for a hundred orbits or so. 

It is tempting to suggest that the observed shape of the disk rim of HD 142527 today 
will also be brief compared to its age, with features that are unlikely to be found in other 
systems that are also undergoing dynamical clearing. In this scenario, snapshot imaging 
of transition disks are likely to look different in their fine details. However, there should 
be common features like off-centering of the central star, non axi-symmetry of the disk 
rim, and presence of spiral a rms. Current imag ing campaigns are indeed reveali ng these 



structures, as in 



^Hashimoto et al. 



D 100546 (IGradv et al. 



200l|); HD 135344B flMuto et all 120121 ); AB Aur 



20111 ). Interstingly, the list of known transition disks is rapdily increasing 



with improving observational capabilities both in dir ect imaging; at optica l/near-IR (see 



e.g., the citations above) and in the sub-mm regime ( [Andrews et al 



201l|). We should be 



able to verify these ideas in the near future, specifically that the observed morphological 
variety of transition disks derives from a common dynamical history in this evolutionnary 
stage when there is large-scale feedback between the planet formation process and its parent 
disk. 

A key to test this transient scenario for HD 142527 and transition disks in general 
would be to detect the gap-crossing planetary accretion flow and link it to the outer-disk 
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asymmetries. Meanwhile, a quantitative comparison with hydrodynamic models requires 
radiative transfer calculations on the model density fields. 



While this article was under review 



Biller et al. 



( 120121 ) reported on a tentative detection 



of a low-mass stellar companion to HD 142527, at ~13 AU, so abutting on the inner 
disk. This companion is inferred through parametrised- mo deling of visibility data. If 
real, this stellar companion will strongly impact on the inner disk, and its inclusion in 
simulation may bring closer agreement with the observed inner disk properties. This would 
be particularly interesting given that our one-planet scenario currently fails to reproduce 
the properties of the inner disk. It is currently too large in radius and mass, compared with 
observations. Including a second body in the simulation w ould produce a larger gap with 



a smaller inner disk (e.g, 



Dodson-Robinson fc Salyk 



201 lh . However, HD 142527B alone 



probably cannot explain the extent of the cavity. A single body can not clear a gap wider 



Dodson-Robinson fe Salyk 



201 lh . given by Rmi, where 



than ~5 times its tidal radius (e.g. 
-Rhiii — fl(^p/3M*) 1//3 is known as the Hill r a dius. Mass and semi-major axis estimates for 



the companion star detected by 



Biller et al. 



J2012|) give a Hill radius R m - 3 - 5 AU. 
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1 0-11 0-1 




1 0-11 0-1 



Fig. 1. — Infrared adaptive-optics imaging of HD 142527 and its protoplanetary disks, from 
ll-Jun-2011. In x and y we indicate offset J2000 RA & Dec, in arcsec (so North is up and 
East is to the left). From left to right we show reduced images obtained in L p , H 2 (1-0)S(1), 
and K s , and an RGB combination following the color codes indicated in each image. The 
linear gray scales are normalised over the entire range of intensity values outside a radius of 
0.45 arcsec, containing the core of the diffracted stellar light, and abutting on the halo of the 
stellar glare. Inside this region the images are masked, except for the H2 image, where the 
stellar PSF has been attenuated with azimuthal averages. A roughly elliptical outer ring is 
clearly seen. The star is offset from the center of the cavity by 0.17±0.02 arcsec. 
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1 1 

Fig. 2. — Annotations on the H 2 (1-0)S(1) image also show on Fig. [TJ In red we indicate the 
spiral arms we could identify, with red identification numbers. The green numbers indicate 
the positions of the two intensity nulls, or gaps, along the ring. The blue box surrounds the 
position of the knot. 



-14- 




0.2 -0.2 

Fig. 3. — Upper limits on the K s magnitude of HD 142527b. Axis labels are as for Fig. [TJ In 
a) (left panel) we show the collapsed SINFONI datacube, as produced by the pipeline, but 
divided by an azimuthal average chosen to highlight the PSF pattern. In b) (right panel) we 
show the collapsed image of the spectrally-deconvolved and PSF-subtracted data cube, with 
a linear intensity scale streching from -0.75 to 0.57, in counts, where the maximum number 
of counts in the collapsed cube (shown in a)) is 548. The root-mean-square dispersion of the 
counts in an annulus from 0.2 to 0.3 arcsec is 9.6 10 -2 , or a 3 a lower limit to the flux ratio 
of 1905. Using a separate PSF standard, this residual rms value is 0.19, giving a 3 cr lower 
limit to the flux ratio of 966. 
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Fig. 4. — FARGO simulation of a 10 Mj up planet embedded in a 0.1 M© disk. The spatial 
axes are in AU. The color scale represents the surface density in a logarithmic ramp in 
FARGO code units in which 10~ 5 equals 88 g cm -2 . Notice the stellar shift relative to the 
position of the cavity centroid of ~16 AU. 
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